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Abstract 

Two different gauge potential methods are engaged to calculate explicitly the spin Hall conductivity 
in graphene. The graphene Hamiltonian with spin-orbit interaction is expressed in terms of kinematic 
momenta by introducing a gauge potential. A formulation of the spin Hall conductivity is established 
by requiring that the time evolution of this kinematic momentum vector vanishes. We then calculated 
the conductivity employing the Berry gauge fields. We show that both of the gauge fields can be 
deduced from the pure gauge field arising from the Foldy-Wouthuysen transformations. 
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1 Introduction 



For low energies and at the Dirac points where conduction and valance bands meet, graphene is 
described by some copies of the two-dimensional massless, free Dirac Hamiltonian after substituting 
the Fermi velocity Vp for the velocity of light c [1]. Mass gaps can be generated by taking into account 
spin-orbit interactions. One of the interesting features of this Dirac-like theory is the accomplishment 
of the spin Hall effect in graphene [2j. It is based on the observation that within this Dirac-like theory 
the quantum Hall effect occurs in the absence of external magnetic fields [3j. This discovery of Kane 
and Mele (KM) [2] had far reaching consequences like a new phenomena called topological insulators 
(see [4] for a complete list of references). We would like to consider the Dirac-like Hamiltonian of 
KM and calculate the spin Hall conductivity in graphene directly using quantum mechanical and 
semiclassical methods. 

We will present two different but somehow related approaches. The first method relies on the 
observation that the Hamiltonian proposed by KM can be rephrased by defining a non-Abelian gauge 
potential. We will deal mainly with the spin-orbit interaction which preserves spin. When the Hamil- 
tonian is expressed in terms of a gauge potential one can introduce the related kinematic momentum 
vector and define its time evolution which can be conceived as the quantum mechanical analogue of the 
force. In the case of the Hall effect this is the Lorentz force. One may provide a method of calculating 
the Hall conductivity through the configurations obtained by requiring that the Lorentz force vanishes 
which can be generalized to the integer quantum Hall effect. We will show that acquainted with the 
gauge potential it is possible to derive explicitly the spin Hall conductivity in graphene imitating this 
formulation. It is in accord with the result of KM. 

The Foldy-Wouthuysen transformations of the two-dimensional massive Dirac Hamiltonian can 
be employed to define a non-Abelian gauge fieldjS]. Because of being pure gauge field it leads to 
a vanishing field strength. However, in the adiabatic limit it generates the Abelian Berry gauge 
potential[6j. Examining the semiclassical equations of motion of electrons in the presence of Berry 
gauge fields, one can show that the Hall conductivity is related to the Berry curvature. We generalize 
this formulation to the Dirac-like Hamiltonian of KM and show explicitly that spin Hall conductivity 
can be written in terms of the Berry phases. On the other hand we will show that the gauge potential 
of the KM Hamiltonian which is furnished in the first method can be deduced from this non-Abelian 
gauge field in the limit of vanishing momenta. 

The method of calculating the Hall conductivity through the configurations defined by the require- 
ment that the time evolution of the canonical momentum vector vanishes is presented in Section [2j 
In Section [3] we first express the KM Hamiltonian in terms of a vector field. We then imitate the 
formulation of Section [2] to derive the spin Hall effect conductivity in graphene by taking into account 
only the spin conserving spin-orbit interaction. Section U] is devoted to the calculation of the spin Hall 
effect in graphene in terms of the Berry phases. In the last section we discuss the relation between 
these two different gauge fields and comment on whether the other possible interaction terms can be 
treated similarly. 

2 Integer quantum Hall effect 

To illustrate the formalism which we will adopt to calculate the spin Hall conductivity in graphene 
we first would like to discuss the integer quantum Hall effect. Hall effect is a phenomena which 
occurs when electrons are constrained to move in a two-dimensional sample which is subjected to 
the uniform, perpendicular magnetic field B. Because of the Lorentz force the charge carriers will 
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move to one edge of the sample creating a potential difference between the two edges. By applying 
the in-plane, constant electric field E one can balance this potential difference thus the electrons will 
move through the sample without deflection. Considering the velocity of the electrons in this system 
one can find the Hall current and derive the Hall conductivity [7]. It is possible to carry this approach 
into quantum mechanics and obtain the integer quantum Hall effect. To this aim let us introduce 
the kinematic momenta Hi = pi — a, associated with the electromagnetic vector potential components 
ai = {eB /2c)eijrj] i,j = 1,2, and consider the Hamiltonian 

Hs = ^-eE-r. (2.1) 

The commutation relations obeyed by the kinematic momentum and position operators are 

[ri,Tij\ = ihdij, [ri,rj] = 0, [TTi,7rj] = ihfij, 
where fij = daj/dri — dai/drj is the electromagnetic field strength with the nonvanishing components 

/i2 = -/21 = -B. 

c 

The Heisenberg equations of motion following from the Hamiltonian (12. ip are 

Tj = ^, 2.2 

m 

TTi = eEi + fijr-j. (2.3) 

Obviously, the eigenvalues of these operators correspond to the velocity of charge carriers and the 
electromagnetic force acting on them, respectively. Thus quantum mechanically demanding that the 
electrons do not be deflected by the electromagnetic force, is equivalent to consider the states ipnif) 
satisfying 

^V'H(r) = 0. (2.4) 

One can observe that nontrivial ipnir) exists. This ensures that we can consider the subspace where 
()2.4p is satisfied and solve them to obtain the velocities of charge carriers which are not deflected as 

c 

Vi = -^^ijEj. 

Although we do not indicate explicitly these are the values of the "velocity operators" (j2.2p in the 
subspace where the conditions ()2.4p are fulfilled. Plugging these into the definition of the electrical 
current 

j = —CKV 

where k is the concentration of charge carriers, yields 

ji = O'H^ijEj. 

Thus the Hall conductivity is obtained as 

cr. = (2.5) 
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We should also specify k. The eigenstates of the Hamiltonian p.ip for E = 0, which describes 
the electrons moving on the plane in the presence of the perpendicular, constant magnetic field B are 
known as Landau levels. They correspond to the eigenvalues 

En = hwc{n + ^); n = 1, 2, ■ ■ ■ , 

where = The density of states of each Landau level can be calculated to be 

m 



Pl{E) 



For the two-dimensional free electrons one finds the same density of states. One can attain the 
concentration of charge carriers between the the ground state Eq and the Fermi level E^^ as 

^%,{E)dE = ^{E^^^-Eo). (2.6) 

For Landau levels Eq = ^ and Ep^^ = fiiOc{N + |) with N being the index of the highest filled 
Landau level, so that (12. 6p yields 



Inserting ()2.7p into ()2.5p one obtains 
which characterizes the integer quantum Hall effect. 
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3 A gauge field formulation of graphene 

Graphene is the one atom thick material composed of carbon atoms which are situated at the corners 
of hexagons arranged as a honeycomb lattice. It is constructed as a superposition of two triangular 
sublattices which are designated as A and B. Its conduction and valance bands touch at two inequiv- 
alent Dirac points K and K' . Around these points, for low energies, the particles are described by the 
free, massless, two-dimensional Dirac-like Hamiltonian [1] 

HQ = VpCx-ip, (3.1) 

where Vp is the effective velocity of electrons and 

oi. = (axT-^,ay) . (3.2) 

^x,y,z are the Pauli matrices acting on the states of the sublattices A and i?, in the representation 
where Oz = diag(l,— 1). The other Pauli matrix = diag(l,— 1) labels states at the Dirac points 
K and K' . One can also introduce the Pauli matrices Sx,y,z to identify the spin of electrons with 
Sz = diag(l, —1). We do not explicitly state neither the direct product between these spaces nor the 
unit matrices. In the presence of the uniform magnetic field B in z direction the spin-orbit interactions 
proposed in [2] yield the Hamiltonian 

eB 

H° = Ho + Aso<yzTzSz + XnicTxTzSy - GySx) ^(ycr^' " XGyTz)- (3.3) 
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To imitate the formalism of Section [2] we need to express the interaction Hamiltonian ()3.3p by means 
of gauge fields. Obviously, the magnetic field term can be written in terms of an electromagnetic gauge 
field. However, the fully fledged interaction Hamiltonian can be written as 

H° = vpcx ■ 11° 

by substituting p in the free Hamiltonian (j3.ip with the kinematic momentum vector 

by introducing the following gauge potential components 

A- = ^^'^y^^-^^y + ^/^y^ (3-4) 



AG -^so , A^j eB 
A = -I- — axTzSz H Sx - - — T^x. (3.5) 

" 2vp Vp 2vp 



The last terms are the ordinary electromagnetic gauge fields in the symmetric gauge and the second 
terms can be read directly from ()3.3p . The first terms are unusual and they arise through the algebra 



of the Pauli matrices: axfJy = icFz- 

The operators H^ and rj satisfy the commutation relations 

h, Hf] = ih5i,, [ri,rj] = 0, [Hf , Hj] = ihet^F^y. (3.6) 

The field strength is defined as 

dA'^ dA? i 

Thus, for the non-Abelian gauge field (j3.4p . (j3.5p . it takes the form F^- = €ijF^ with 



dr, dr, h^^''^^^- 



eB iAcoA„, , 2A? 



Fx% = -r. - ^a.r. + -^{aysy + axT.Sx) + (3.7) 



3.1 Spin Hall effect in graphene 

To discuss spin Hall effect we require that only the the spin-orbit term conserving the third component 
of spin is nonvanishing in (j3.3p . Hence, we deal with the Hamiltonian 

H = Ho + AsofT,T,s„ (3.8) 

which can written as 

H = Vpot ■ n, 

through the kinematic momenta Tl = p — A, where the gauge field A is obtained from (j3.4p - (j3.5p by 
setting B = 0, = : 

Ax = i— o'ySz, Ay = o'xTzSz- (3-9) 



This is a vector potential taking values in the SU{2) group generated by (yx,y,z- The curvature corre- 
sponding to (13. 9p is 



'-^so 



Fxy- ^ \Ax,Ay] - ^^nCTzT, 
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We are now equipped with the whole machinery needed to discuss the spin Hall effect imitating the 
formulation of the Hall conductivity presented in Section 2. Considering the electrons in the external 
constant electric field E described by 

HE = H-eE-r, (3.10) 

we will derive the Heisenberg equations of motion and obtain the spin current by requiring that the 
time evolution of the kinematic momentum vector vanishes. The "velocity" and "force" operators 
obtained through the Heisenberg equations of motion r = j^[HE,r] and 11 = ^[//b,!!] are 



ri = {x,y) = {vpa^T;,,Vpay), (3.11) 
iA 

ili = -^^CTzTzSz^i - ^^eijrja;,Tz + eEi. (3.12) 
Then, we demand that the "force" vanishes: 

Uqf{r) = 0. (3.13) 

One can easily show that for all of the components of the spinor wave function ^'(r) the same differ- 
ential equation follows from (|3.13p : 



f^W = /W(l, 1,1, 1,1, 1,1,1), (3.14) 

f{r) = 0, (3.15) 



A2 A"^ 
-j^P. + {eE,) + ^ 



where i = 1, 2 is not summed over. In the rest of this section we deal with operators taking values in 
the subspace spanned with the spinors (j3.14p and (|3.15|) . although we will not explicitly write. In fact 
in this subspace we can now solve ()3.13p to obtain the velocities 



2ivl 



eijUjSz + rHi, (3.16) 



^so 

where we separated the electric field dependent part as 

2ehvleijEj 

rm = .n CTzTz- (3.17) 

^so 

We would like to calculate the spin Hall conductivity, therefore the relevant parts in the solution 
(j3.16p are the terms proportional to the electrical field components (j3.17p . Let us label the velocity 
of spin up carriers by r^- and the velocity of spin down carriers by rjj-. We can further determine the 
velocities of particles in the K and K' valleys corresponding to 1 and —1 eigenvalues of as 



2ehvleijEj _ 2efv4ei^ 

'-^so '-^so 

4+ 2ehvleijEj 4_ 2ehvleijEj 

''Hi = .2 = A2 (^-19) 

The Hall currents of the spin-up and spin-down electrons are defined as 

jl = n^+rl+ + n^-rl-, (3.20) 

jj^ = n^+ri^ +n^-ri~, (3.21) 
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where n^~^,n^ ,n'^~^,n'^ indicate concentrations of the related carriers, The Hall currents (|3.20|) . 
()3.2ip can be employed to define the spin Hall current as 



Jh 



H Jh) 



(3.22) 



We need to determine the concentrations n^~^,n^ ,n-^'^,n'^ . This will be elaborated inspecting the 
corresponding Hamiltonians. In fact, there are four different two-dimensional Hamiltonians stemming 
from dlSl): 

/ \ 

H^- 

V H^- J 

These two-dimensional Hamiltonians corresponding to the ti i spin and the K, K' valley are 



H 



H^^ = Vp{axPx + CTyPy) + AsoCTz, 
= Vp{axPx + CTyPy) - AsoCTz, 



= Vf{-(TxPx + CTyPy) - AsoCTz 

= Vp{-axPx + CTyPy) + AsoCTz 



(3.23) 
(3.24) 



The effect of the spin-orbit term is to create a gap in the energy band structure of the Hamiltonians. 
In terms of the eigenvalues of the momenta hk, (j3.23p and (|3.24p yield the same energy distribution 



(3.25) 



corresponding to particle and antiparticle (hole) states. We choose A50 > and let the Fermi energy 
of graphene be in the gap by setting A; = 0. Identifying the concentration of particles and antiparticles 
(holes) by Up and n^, we obtain 



n 



t+ 



Up 

Ha 



n 







Ur, 



n' 



rv 



Ua 
n„ 



Up 

Ua 



(3.26) 
(3.27) 



The derivations of (|3.26p and (|3.27p are elaborated in Appendix A. Inserting the carrier concentrations 
(|3.26p and (j3.27p into the definition of the spin Hall current (I3.22p one obtains 



J HI 



(3.28) 



where n = Up — ria- Now, n can be calculated using the density of states corresponding to the energy 
distribution ()3.25p found as 



Pe{E) 



(Pk' 



\E\ 



The edges of band gap lie at Aso and —Ago, hence when we consider positive Fermi energy which 
lies in the gap, the number of states is derived as 



n 



Aso 



PE{E)dE 



\Aso\Aso 



(3.29) 



7 



Plugging (Km into (pTM]) yields 



3 Hi 




and we obtain the spin Hall conductivity as 



This is in accord with the results of O [2] . 

4 Spin Hall conductivity as Berry phase 

Topological nature of the Hall effect is well exhibited in terms of Berry phases. The semiclassical 
equations of motion are altered drastically in the presence of Berry gauge fields (see and the 
references therein) . They yield an anomalous velocity term for electrons which leads to the anomalous 
Hall conductivity. In fact, ignoring spin of electrons the Hall conductivity can be written in terms 
of the Berry curvature Fg on the Fermi surface [9l [10] as (a complete list of references for the Berry 
phase effects in this context can be found in the recent review [llj ) 



Considering the electrons with spin, a generalization of (j4.ip to the spin Hall effect was discussed in 
[T^ . Note that in this section p is not a quantum operator but denotes the classical phase space 
variable. We deal with four different two-dimensional Dirac-like theories (j3.23p . (j3.24p . thus we should 
take into account the contributions arising from each of them separately. We adopt the formulation 
of [5] to derive the Berry gauge fields arising from each one of the two-dimensional Hamiltonians 
(13.23p . (13.24p . Therefore, we should start with giving the unitary Foldy-Wouthuysen transformations 
(C/^"*", J/^"*", U'^~) corresponding to the Dirac-like Hamiltonians (-ff^"*", , H-^'^ , H-^~). We would 
like to present them in the unified notation: 



The unitary Foldy-Wouthuysen transformation U can be engaged to define the gauge field|131 [5] 



Exploring the Dirac-like Hamiltonians (|3.23p , (|3.24p , we introduce the following Foldy-Wouthuysen 
transformation 




(4.1) 



U = diag(?7^+, U^-, C/^+, U^-). 



A = diag{A'^-^,A'^-,A^+,A^~) = ihU{p) 



dU^p) 
dp 



( a,H^+ + E 





<j,H^- + E 









cr^i7-^+ + E 




\ 



u = 



1 












a,H^- +E j 



(4.2) 



^2E{E + J\so) 



\ 




Y^f^p^ + A|q. Observe that (|4.2p is defined to 



(4.3) 
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One can study each entry of (14. 2 p as in [5J and show that they lead to the gauge potential 
ih 

= oF2^F I A — \ [vfE{E + Aso)a(TzTzSz + vla.,T-,s^{a ■ p)p + vlE{a. ■ p)a - vlpE] . (4.4) 
Its components can be written explicitly as 

. _ h{VpE{E + Aso)aySz - vljcTyPx - CFxTzPy)SzPx + vlEa-^T-^Py) 

^ " 2E^{E + Aso) ' ^ ^ 

. _ h{-VpE{E + /\so)(yxTzSz - vljcTyPx - CrxTzPy)SzPy - viEazTzPx) . 

' ~ 2E\E + Aso) • ^ ^ 

Because of being a pure gauge potential the field strength of (j4.4p vanishes. However, one can consider 
the adiabatic approximation by projecting on the positive energy states: 

= PAP. (4.7) 
By inspecting the positive eigenvalues of the ()4.3p one can deduce that the projection operator P is 

P = diag(P;, pI, P:^' ) = diag(l, 0, 0, 1, 0, 1, 1, 0). 



It should be noted that when projected on positive energy states only the last terms in (|4.5p and (j4.6p 
make nonvanishing contributions, so that the Abelian Berry gauge field is 

= 2EiE ^Asof -'^^^'- ^'-'^ 

where the unit matrix 1^ in the space is exhibited explicitly. It is worth emphasizing that in A 
negative energy states are present, thus it possesses twice the matrix elements of A^ . The nonvanishing 
component of the Berry curvature is given as 



We propose to generalize (|4.ip to the spin Hall effect in graphene as 

+ Tl)-{J^i+ Ti)] . (4.10) 



g e f F dp 



" 2 J {2TTh) 

1 

into the definition (|4.10p leads to 



Ep denotes the highest energy level occupied in the two-dimensional system. Thus, inserting ()4.9p 



S _ e / - u p I VpLJS.so 



d?p ( 2f?v\A 



2 (2'KhY V E^ 

e Ago 



(4.11) 



2vr • 

This in accord with the calculation of the spin Hall conductivity obtained by employing the Kubo 
formula which is presented in Appendix B. 

/ox 

We let the Fermi energy level of graphene lie in the gap, so that in (14. lip we set Ep — Ago and 
obtain the spin Hall conductivity as 



27r 

This is the value established in [2]. 
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5 Discussions 



The main difference between the gauge fields (j3.9p and (|4.8p is the fact that the latter is acquired 
from (14. 4|) in the adiabatic limit ()4.7p . thus it is within the particle bands of the two-dimensional 
Dirac-like theories, but the former one is non-Abelian which connects the particle and hole bands of 
the two-dimensional Dirac-like theories. However, also the non-Abelian gauge potential (j3.9p can be 
obtained from the gauge field ()4.4p in the vanishing momentum limit, up to a constant , as 



■\p=0- 



This limit corresponds to restrict the Fermi energy of graphene to lie in the gap. In fact, the calculation 
of the spin Hall conductivity presented in Section 13.11 is valid only for the Fermi energy lying in the 
gap. Therefore, although the calculation methods are different the gauge fields (|3.9p and (|4.8p result 
from the same gauge field (j4.4p in two different limits. 

The nonvanishing magnetic field B ^ case can be studied similarly. However if we switch on the 
Rashba term A/j 7^ in (j3.3p . the third component of spin is not conserved and in spite of the fact 
that we can introduce the gauge potential (j3.4p . (j3.5p . the related curvature (j3.7p is not diagonal in the 
spin space. Hence the method of Section [3.11 is not anymore suitable to discuss the spin Hall effect. 



Appendix A 

Introducing the label / = (t+j t~) 4+) i~) corresponding to each of the two-dimensional Dirac-like 
Hamiltonians given in ()3.23p . (I3.24p . we define the number operators as 

= np\p'){/\+na\a'){a'\, (A.l) 



where \p^) and \a^) denote, respectively, the positive energy E = ^ v%T?k'^ -|- A^q and the negative 
energy —E eigenspinors of the related Dirac-like Hamiltonians. For instance, let us consider the 
Hamiltonian for the spin up carriers in the K valley, -ff^"*", whose eigenspinors can be written in the 
chiral basis as 

where cos0 = ^W- and tanc/) = Plugging them into (jA.ip leads to 



Art+rt.^ - ^^4^ + ^^2^e-*<^sin0 \ 

n£z!^ei<^sine !iE+!i£ _ !^ cos 7 ' ^ ' 

When the Fermi level of graphene is in the gap generated by the spin-orbit interaction we should set 
fc = 0, hence cos B = \ and (IXil) yields 
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In order to obtain the concentrations corresponding to the other Hamiltonians, we need to consider 
their eigenspinors. In fact, for the spin up carriers in the K' vahey described with H^^ the eigenspinors 
are 

,t-\ _ f sin(2, 



1^'") = ( -cos( g;-^^ )' (^-^^ 



sin(|)e 



§.\p-i<t> 



(A.6) 



Inserting them into (jA.ip and considering the vanishing momentum hmit leads to 



N^-{k = 0) 



n 







( sin(f) 
V cos(|)e^'^ 








( cos(|) 
V -sm(|)e^' 







where is given in (|3.26|) . Similarly, for the spin down carriers in the K valley, we can show that 
the eigenspinors of the Hamiltonian H^'^ are 

(A.7) 
(A.8) 

Making use of them in (jA.ip and setting = yield the concentration number n^~^ given in (|3.27p : 

N^+{k = 0) = n^+. 

The eigenspinors of the Hamiltonian H^^ corresponding to the spin down carriers in the K' valley are 



= Z^elUs , (A.9) 



\-'-) = ( JS§.. ) • (A.10) 



cos(|je 



One can easily observe that they result in 



N^-{k = 0) 



n 



4- 



where is given in (j3.27p . 



Appendix B 

Kubo formula corresponding to the Hamiltonians (I3.23[) . (13.24p can be written in the notation of 
Appendix A as[l4j 

. eh^ C^'f <fk 2lm[{a'\y\p'){p'\i\a')] 

= — j (2^ ' ^^-^^ 

where x, y are the related velocity operators which can be read from (|3.1ip and 
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For H"^^ the velocity operators are x = Vpa^, y = VpUy and the eigenspinors \p^^) and \a^^) are given 
in (|X2]) and (|X3]) . Employing them in (|RT]) leads top^ 



For the spin up carriers in the K' valley, we set x = —VpUx-, y = VpGy and deal with the eigenspinors 
()A.5p . ()A.6p . We obtain the same conductivity 

The contributions arising from the spin down carriers in the K and K' valleys are also equal but differ 
in sign with the spin up contributions: 

i<r = K)^- = (B.4) 

To obtain the spin Hall conductivity we should take the difference of the spin up and spin down 
contributions as 

^'u = ((^^)^+ + {<y%r) - {{<y%f^ + {<yl)^-)- (b.s) 

Inserting ([Ell), <lB^ and (iRil) into (IB3]) leads to the spin Hall conductivity 

s e Ago 

This is the same with the result obtained in terms of the Berry phase (14. lip . 
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